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ABSTRACT 

Complete modeling of radiative transfer in neutron star atmospheres is in progress, taking into account 
the anisotropy induced by magnetic fields, non-ideal effects and general relativity. As part of our mod- 
eling, we present a novel atomic calculation method producing an extensive atomic data set including 
energy values and oscillator strengths in the so-called Landau regime (B > 4.7 x 10 9 Z 2 G). Conven- 
tional atmosphere models for B = are not applicable to typical field strengths of cooling neutron stars 
(B ~ 10 12 — 10 13 G), since an atom no longer keeps its spherical shape. The elemental composition 
and the configuration of the magnetic field in the atmosphere are presently unknown, so that atomic 
data must be produced for ground and excited states of several ions as a function of magnetic field. To 
accomplish this efficiently, we minimized the iterations in the Hartree equation and treated exchange 
terms and higher Landau states by perturbation methods. This method has the effect of reducing the 
computation time significantly. Inclusion of higher Landau states gives us much more accurate data for 
inner orbitals unlike other methods based on the adiabatic approximation. While existing atomic data 
in the Landau regime are available only for low Z atoms, our method can be used in elements up to 
Fe with sufficient accuracy to be of use for spectroscopic missions such as Chandra, XMM-Newton and 
next-generation X-ray telescopes. 

Subject headings: atomic processes — magnetic fields — stars: neutron 



1. INTRODUCTION 

Since their discovery in 1967 (Hewish et al. 1968), pul- 
sars have attracted the interest of researchers in various 
fields due to their exotic properties and extreme physical 
conditions. Although many of the pulsars were discovered 
in the radio band, multi-wavelength studies have helped 
to elucidate their properties. ROSAT discovered thermal 
components (T ~ 10 5 — 10 6 K) from the so-called cooling 
neutron stars in the soft X-ray band (Becker & Trumper 
1997). Emergent spectra from the photosphere can be 
modified from a pure blackbody spectrum by the presence 
of an atmosphere (Romani 1987). The geometrical thick- 
ness of neutron star atmospheres is small (~ 0.1 — 100 
cm) due to the strong gravitational field (Pavlov et al. 
1995). Nevertheless, it can be optically thick for repro- 
cessing thermal photons due to the high density. 

With great improvements in the sensitivity and the res- 
olution of recent X-ray telescopes, X-ray spectroscopy can 
directly provide information on the surface of neutron 
stars. There are two possible implications that the at- 
mosphere may have on the X-ray spectrum. Precise mea- 
surement of surface temperature taking into account the 
atmospheric effects is crucial for determining the cooling 
curves of neutron stars, which constrains various equations 
of state (Tsuruta 1995; Page 1998; Yakovlev et al. 1999). 
Secondly, atomic (line or edge) absorption features along 
with cyclotron lines enable us to diagnose properties of the 
atmosphere. 

Magnetic field (B) can be measured directly by the ob- 
servation of cyclotron lines (Trumper et al. 1978; Wheaton 
et al. 1979), or indirectly by the spindown rate of pul- 
sars assuming a dipole magnetic configuration (Ostriker & 
Gunn 1969). Most of the observed pulsars are believed to 



possess very strong magnetic fields (B ~ 10 11 — 10 13 G), 
except the millisecond pulsars whose magnetic fields are 
relatively low (B ~ 10 8 - 10 10 G) (Camilo et al. 1994). 
Surface elements (Z) are not known at present. Strong 
gravity stratifies the surface layers so that the lightest 
element resides at the top of the atmosphere (Alcock & 
Illarionov 1980). However, any element between H and 
Fe is feasible. Accretion from the interstellar medium 
could accumulate a hydrogen layer, while an Fe atmo- 
sphere is favored from the study of the mass-cut in su- 
pernova explosions (Woosley & Weaver 1986). Lighter el- 
ements are formed by the spallation of heavy nuclei by 
particles (Miller 1990). On the other hand, pycnonuclear 
reactions could produce heavier nuclei from light nuclei 
(Lai & Salpeter 1997; Salpeter 1998). The fallback from 
the outer shells of supernova remnants could be of any 
element (Chevalier 1996). The surface density (p) is also 
uncertain in the range of p ~ 0.1 — 10 gem - (Pavlov et al. 
1995). The state of matter on the surface could be liquid, 
solid or in the form of molecular chains along the field 
(Abrahams & Shapiro 1991; Lai et al. 1992; Demeur et al. 
1994; Ortiz et al. 1995; Lai & Salpeter 1996, 1997). In high 
density plasma, atoms are strongly-coupled, inducing non- 
ideal effects, and electrons may be partially degenerate. 

2. atmosphere models 

There have been several attempts at modeling magne- 
tized neutron star atmospheres. For weakly-magnetized 
pulsars (B < 10 9 G), atmosphere models for B = are 
applicable (Romani 1987; Rajagopal & Romani 1996; Za- 
vlin et al. 1996). These non-magnetic atmosphere mod- 
els have been applied to several neutron stars (Zavlin & 
Pavlov 1998; Pons et al. 2000). However, for the strongly- 
magnetized atmospheres (B > 10 9 G), atomic models re- 
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quire a new scheme quite different from the conventional 
atomic calculations for the B = case. For instance, 
the ionization threshold of a hydrogen atom increases to 
160 cV at B = 10 12 G. Also, in the presence of a strong 
magnetic field photon opacities are highly dependent on 
the direction of the magnetic field and the polarization 
of photons (Pavlov et al. 1995). Thermodynamic proper- 
ties of hydrogen atmospheres have been intensively inves- 
tigated by Lai & Salpeter (1995), Steinberg et al. (1998) 
and Potekhin et al. (1999). They have computed the equa- 
tion of state and the degree of ionization assuming LTE 
and taking into account coupling effects such as the mo- 
tional Stark field and the pressure ionization peculiar to a 
strongly-magnetized dense plasma. The hydrogen atmo- 
sphere model of Potekhin et al. (1999) is constructed based 
on the accurate atomic database obtained by Potekhin 
(1998). Recently, Ozel (2001) and Ho & Lai (2001) stud- 
ied spectral signatures from atmospheres of strongly mag- 
netized neutron stars, utilizing improvements in radiative 
transfer theory. However, they assumed fully-ionized hy- 
drogen atmospheres. Therefore their results are limited 
to rather high temperature. Similarly, Zane et al. (2000) 
studied fully-ionized hydrogen atmospheres of accreting 
neutron stars. Fully-ionized hydrogen atmosphere mod- 
els have been applied to data on isolated neutron stars 
obtained by XMM-Newton and Chandra (Paerels et al. 
2001; Pavlov et al. 2001). In contrast, radiative transfer 
in partially-ionized hydrogen atmosphere with full inclu- 
sion of bound-bound and bound-free opacities remains to 
be solved in comparison with recent X-ray data. 

On the other hand, atmosphere models of high Z 
atoms are still immature due to their increased complex- 
ity. Miller (1992) was the first to construct model atmo- 
spheres for Z = 2 — 7. However, the opacities include 
only bound-free cross sections with averaging of the po- 
larization modes. Fe atmospheres were investigated by 
Rajagopal et al. (1997); the model is rather crude and 
the energy values and oscillator strengths have as much 
as 10% and a factor of two uncertainties, respectively. 
Both of the atmosphere models for Z > 1 elements are 
based on atomic calculation by the Hartree-Fock method 
of Neuhauser et al. (1987). 

As the spectroscopy capability of X-ray telescopes im- 
proves, model atmospheres of strongly- magnetized neu- 
tron stars should comparably improve. Our efforts are 
focused on constructing atmosphere models for Z > 1 
atoms for comparison with existing hydrogen atmosphere 
models and new X-ray observations. As part of this ef- 
fort, we present a novel atomic calculation appropriate for 
model atmospheres of strongly-magnetized neutron stars. 
A brief review of atomic structure in strong magnetic fields 
is given in §3, then various atomic calculations are de- 
scribed in §4 in comparison with our formalism in §5. Sub- 
sequently, the validity of our model and relevant physical 
processes are discussed in §6. Some numerical improve- 
ments are mentioned in §7. Finally, we present our results 
in comparison with previous atomic models in §8. 

3. ATOMIC STRUCTURE IN STRONG MAGNETIC FIELD 

When a strong magnetic field is present, the spherical 
symmetry of an atom is broken and the atom is stretched 
along the field. The Landau regime is defined as that 



where the magnetic field effects exceed Coulomb field ef- 
fects, i.e. fkdB > Ze 2 /r (Hujb = ll.57.Bi2 keV = electron 
cyclotron energy, B = 10 12 Bi2 G). This is translated to 
j3 z > 1, where j3 z = B/Z 2 B , f3 = B/B a and the refer- 
ence field Bq = 4.701 x 10 9 G, as a good measure for the 
relative strength of magnetic field to the nuclear Coulomb 
field. The dominance of the magnetic field requires us 
to investigate the atomic structure exploiting cylindrical 
symmetry. In cylindrical coordinates (p, <fi, z) bound elec- 
trons in an atom are well described by the four quantum 
numbers (n,m,v, s z ) denoting Landau number, magnetic 
quantum number, longitudinal quantum number and com- 
ponent of electron spin along the field respectively, v is the 
node of the longitudinal wavefunction and it defines the 
parity along the field as (— X) v . v = states are located 
close to the nucleus (tightly-bound orbitals), while v > 
states are far from the nucleus (loosely-bound orbitals). 

In very strong magnetic fields (f3z 3> 1) it is a fairly good 
approximation to fix n = (adiabatic approximation) and 
m < 0, s z = —1/2 (full-spin-polarization, FSP hereafter). 
Various atomic models have adopted both assumptions 
(i.e. restrict the two quantum numbers n and s z ) because 
they simplify calculations. Accordingly, bound electrons 
are characterized by a set of two quantum numbers (m, v) 
(m indicates its absolute value, hereafter). However, as 
Pz approaches 1 the occupation in n > Landau levels is 
not negligible (configuration mixing), so the adiabatic ap- 
proximation does not provide accurate atomic data (Jones 
et al. 1999b). With further decrease in the magnetic field, 
spin-flip transitions take place, i.e. the FSP structure is no 
longer a ground state since an electron with the opposite 
spin lowers the total energy (Ivanov & Schmclchcr 2000). 

The intermediate field regime is defined as that at which 
the electric field and magnetic field becomes comparable. 
Bound electrons feel differing electric fields depending on 
their distance from the nucleus and the screening of nu- 
clear charge by inner electrons. A critical field strength 
B c which defines the intermediate regime is determined 
by equating the Coulomb force and the Lorentz force. 
Assuming that the distance of the (m, 0) orbital to the 
nucleus is ~ (2m + l) 1 / 2 /? where p is the cyclotron ra- 
dius p = 2.566 x 10 _10 i? 12 1 ^ 2 cm and Z e ff is the effective 
charge, B c is given by (Miller 1990), 



A choice of either pure spherical or cylindrical ba- 
sis functions is not sufficient for accurate atomic struc- 
ture calculations in the intermediate field regime. Non- 
separability of p and z in the nuclear Coulomb term pre- 
vents an analytical solution of the Schrodinger equation 
even for hydrogenic atoms. For Z > 1 atoms, electron- 
electron interactions add more complication to the prob- 
lem. 

4. ATOMIC CALCULATIONS IN THE LANDAU REGIME 

Various atomic models have been proposed in the Lan- 
dau regime. First of all no experimental data is available at 
high magnetic fields B > 10 7 G. Therefore atomic data are 
highly model-dependent. Hydrogen has been investigated 
by various authors for the last 20 years (see Ruder et al. 
(1994), Lai (2000) and references therein). Very accurate 
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energy values and oscillator strengths were obtained at ar- 
bitrary field strength with full inclusion of various effects 
such as relativistic correction (Chen & Goldman 1992) and 
motional Stark effect (Potckhin 1998). 

For Z > 1 atoms, the number of studies is signif- 
icantly smaller. Atomic calculations for multi-electron 
atoms are often classified as ab initio or self-consistent 
methods. Ab initio methods include the restricted vari- 
ational method (Chen et al. 1974; Flowers et al. 1977; 
Glasser & Kaplan 1975; Glasser 1975; Miiller 1984) and 
statistical methods. The restricted variational method pa- 
rameterizes the wavefunctions in an ab initio way, so the 
accuracy usually remains poor compared with other meth- 
ods. Statistical approaches such as the Thomas-Fermi 
method (TF hereafter) or Thomas-Fermi-Dirac method 
(TFD hereafter) (Skjervold & Ostgaard 1984; Lieb et al. 
1992; Rognvaldsson et al. 1993; Thorolfsson et al. 1998) 
are applicable to an atom with many electrons, usually for 
Z > 10. 

In general self-consistent methods are more accurate 
than ab initio methods, but they are computationally 
intensive since they entail iterative calculations. The 
one-dimensional Hartree-Fock method of Neuhauscr et al. 
(1987) (1DHF hereafter) was the first self-consistent 
method for multi-electron atoms in the Landau regime. 
Since the adiabatic approximation is assumed in their cal- 
culations, the Hartree-Fock equation reduces to a one- 
dimensional equation along the magnetic field. They pro- 
vided energy values of ground states for Z = 1 — 18 and 
26. Miller & Ncuhauser (1991) obtained energy values 
and oscillator strengths for Z < 14 based on Neuhauser's 
atomic code. As mentioned by Miller & Neuhauser (1991), 
the inclusion of n > Landau levels to the Hartree-Fock 
equation is extremely difficult. 

In recent years several sophisticated atomic models were 
proposed and studied for low Z atoms (Jones et al. 1999b; 
Ivanov & Schmelcher 1998, 1999, 2000). One of them is the 
two-dimensional Hartree-Fock method (2DHF hereafter), 
which unrestricted wavefunctions both in the transverse 
and longitudinal direction (Ivanov & Schmelcher 2000) . It 
unrestricted the electron spin as well. Therefore, 2DHF 
is in principle applicable at any field strength, providing 
very accurate results through improved numerical accu- 
racy. Jones et al. (1999b), who also adopted a similar un- 
restricted Hartree-Fock method with even-tempered gaus- 
sian basis functions (UHF hereafter), pointed out that the 
adiabatic approximation provides accurate atomic data 
only at (iz > 50. Since typical magnetic fields of neu- 
tron stars can be at (3z < 50 for high Z surface elements, 
energy values obtained by 1DHF could be severely inaccu- 
rate in the soft X-ray band where high energy resolution 
spectroscopy is now available. On the other hand, atomic 
data from 2DHF and UHF were computed only for en- 
ergy levels of ground states and several low-lying states of 
low Z atoms (Z < 6 for UHF, Z < 10 for 2DHF) due to 
computation time limits. However, spectroscopy requires 
energy levels and oscillator strengths for various excited 
states as well. 

The density-functional method (DF hereafter) was ap- 
plied to multi-electron atoms in the Landau regime (Jones 
1985a,b, 1986; Kossl et al. 1988; Relovsky & Ruder 1996; 
Johnson & Yngvason 1996; Holas & March 1997). It pa- 



rameterizes the relevant terms in the Hamiltonian (includ- 
ing the exchange term and electron correlation term) as a 
function of electron density in the Kohn-Sham equation. 
However, results are highly dependent on the choice of 
the exchange-correlation density functional (Jones 1985a). 
Ncuhauser et al. (1987) mentioned that the density func- 
tionals of Jones (1985a, 1986) and Kossl et al. (1988) are 
not quite correct in high magnetic fields. In addition, DF 
becomes inaccurate for an atom with a few electrons. 

Other atomic calculations with accuracy ~ 10 -4 such 
as the Quantum Monte-Carlo (QMC) method or the Con- 
figuration Interaction (CI) method have been proposed 
but have been implemented only for two-electron systems 
(Jones et al. 1997; Scrinzi 1998; Becken et al. 1999; Becken 
& Schmelcher 2000). CI method has been extended to high 
angular momentum states of helium (Becken & Schmelcher 
2001). Nevertheless, it required 9 months of computing 
relevant matrix elements to produce an exhaustive atomic 
database. Recently, spectral features from the white dwarf 
GD229 have been successfully explained by helium in a 
magnetic field of fewxlO 8 G (Jordan et al. 1998; Jones 
et al. 1999a). 

Given the lack of experimental data, the self-consistent 
method is required to generate atomic data sufficiently 
accurate for refined X-ray spectroscopy. An extensive and 
accurate atomic data set for various atoms and ion states 
from Z = 1 to 26 in the Landau regime is necessary for 
the X-ray spectroscopy of cooling neutron stars. This con- 
trasts with recent work which has concentrated on low Z 
atoms spanning a wide range of magnetic field, achieved 
with 2DHF and UHF. 

In the next section, after a brief description of our for- 
malism, we describe our novel approach, which we call 
the multi-configurational, perturbative, hybrid, Hartree, 
Hartree-Fock method (MCPH 3 ). The technique allows us 
to achieve an exhaustive list of energy values and oscil- 
lator strengths at various magnetic fields with very fast 
computation times on small computers. 

5. FORMALISM 

5.1. Schrodinger equation preliminaries 

The non-relativistic Hamiltonian in a uniform magnetic 
field (assuming infinite nuclear mass) H is divided into a 
0th order Hamiltonian and a perturbation We 
set the electron g-factor to 2. 




= JT(°)+ff(i). 
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The vector potential of the uniform magnetic field A = 
i_B x f is chosen so that the magnetic field is along the z 

coordinate B = V x A = Bz. Hereafter, all the lengths are 
measured in units of p. V nuc (fi) is the nuclear Coulomb 
potential. V e ff (fj) is the effective potential for electron i 
and represents the mean Coulomb potential from the other 
electrons via the Hartree potential. 



Veff(r i ) = jYl 



3 \zt zt I 



consists of single-orbital Hamiltonians : 



(4) 



(5) 



For an orbital (m, v) (hereafter we omit the index i), the 
Schrodinger equation in the Hartree approximation be- 
comes, 

(6) 



5.2. Multiconfigurational, Perturbative, Hybrid, Hartree, 
Hartree-Fock: The Approach 

Given the formalism above we can now outline the 
MCPH 3 approach and describe its advantage. The 
MCPH 3 method starts with equation (6), which is the 
single-orbital Hartree equation for the unperturbed Hamil- 
tonian H^\ We note that (1) (3 Z > 1 and (2) the 
Hartree potential is orbital independent. The former im- 
plies the occupation of n > Landau levels is small 
compared to the n = state. Therefore equation (6) 
is amenable to a solution in which we break the single- 
orbital Hamiltonian (h^) into a Oth order term (M ^) 
and a 1st order term (h^). The m ' term will rep- 
resent a purely singlc-configurational calculation for the 
single-orbital wavefunction. The orbital-independence of 
the Hartree-potential will ensure that we can restrict the 
iterative calculation of the Hartree wavefunction to the 
n = singlc-configurational (one dimensional) calcula- 
tion; once we have determined the longitudinal part of 
this wavefunction (f mv ) iteratively, it will be immediately 
applicable to the construction of the multi-configurational 
Hartree single-orbital wavefunctions without further need 
for iterative calculations. This represents an enormous 
computational simplification over the schemes employed 
for hydrogen atoms (Rosner et al. 1984; Forster et al. 1984) 
in which the longitudinal wavefunction for each Landau 
level in the multi-configurational solution must be calcu- 
lated separately. 

With a set of eigenfunctions {xnmu) for the Oth order 
singlc-configurational Hamiltonian h^ determined we can 
solve for the multi-configurational single-orbital wavefunc- 
tions {(j) m v) and energy eigenvalues (eml) using perturba- 
tion theory with . We construct the total wavefunction 
ty(r) for H^ by forming a completely antisymmetrized 
spatial wavefunction from the 4> mv (as required by the 
FSP assumption). This total wavefunction \&(r) can then 
be used along with H^ of equation (3) to solve for the 
energy of the total Hamiltonian using simple 1st order per- 
turbation theory. 

This approach is, in effect, a double application of per- 
turbation theory. In the first application perturbation 



theory is used on the single-orbital Hamiltonian to han- 
dle the higher Landau levels and solve for the multi- 
configurational, single-orbital energies and wavefunctions 
(with only 1 iterative calculation required). In the sec- 
ond application 1st order perturbation theory on the total 
wavefunction and total Hamiltonian yields the total en- 
ergy. This same wavefunction is used to calculate oscillator 
strengths. To avoid confusion we call the first application 
perturbation theory of type I and the second application 
perturbation theory of type II. 

This is similar to the Z-expansion method in which the 
whole electron-electron interaction term is treated as a 
perturbation (Hylleraas 1930). However, the Z-expansion 
method is accurate only for highly ionized atoms. Our 
Oth order wavefunction is much more accurate for multi- 
electron atoms due to its self-consistency. As a result, 
we achieved fast convergence (3-5 iterations for less than 
0.1% convergence in total energies) and numerical stabil- 
ity over all the electron configurations we have computed 
for Z = 1 — 26 atoms. The method presented in this arti- 
cle is valid to flz > 1 ; covering the typical magnetic fields 
of cooling neutron stars for the surface elements up to Fe 
(figure 1). 

5.3. The MCPH 3 solution for the multiconfigurational, 
single- orbital Hartree wavefunctions and energies 
( Type I perturbation theory ) 

The single-orbital Hamiltonian h^> in (6) is divided into 
the Oth order and its perturbation as follows. 



where 



(7) 
(8) 



h<® = h B (p, <f>) + h L (z), 
h {1) = {V nuc (f) - V nuc {z)} 

+{Veff(f)-V eff {z)}. (9) 

hB{p,4>) is the transverse Hamiltonian, whose eigenvalue 
and eigenfunction are nhujB and Qnm{p,4>) (see below). 
Due to the separation of variables in h^\ eigenfunctions 
and eigenvalues of the Schrodinger equation for M°) : 

h^ ^ • Xnmi^(^) ^nrav ' Xnmv^J'^) (1^) 

are given by 

(r) = <f> nm (p,<f>)- f mv (z), (11) 

tnmv = nfuV B + eomv (12) 

{Xnmiy} is a complete set of orthogonal eigenfunctions of 
eomv is an energy eigenvalue obtained by the singlc- 
configurational calculation for n = (see below). The 
transverse wavefunction ^ n m(p,<P) is a Landau function 
defined as (Landau & Lifshitz 1965), 



(P \ 

o2 



V/4 



(13) 



satisfying the orthonormal condition 

J pdp j d(j) $* m (/0, 4>) ■ ®n'm'(p,<t>) = 8 nn >5 m m>- (14) 
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5.3.1. Single- configurational calculation for n = 

An equation for h^{z) can be found by multiplying the 
equation (6) by $g m (p, <j)) and integrating over p and <f>. 

h L (z) ■ f mv (z) = e 0mv ■ f mv {z), (15) 

where 

h L {z) = -^^ + V aHe {z) + V. f Az). (16) 

The single-configurational nuclear and effective poten- 
tials arc defined for n = as, 
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V nuc (z) = J pdp J d<t>V nuc (p, z)\$ 0m (p 

V ef f(z) = J pdp J dcf>V eff (p,<P,z)\$0m(p,<P)\ 2 - (17) 



,„,„ with the node v under the 
at z = ±oo with the so- 



We seek an eigenvalue eo 
boundary condition f mv (z) 
lution for f mu (z). The only iterative part to the entire 
solution for H is complete with this calculation toward 
the convergence of eo m v and f mv (z); the one-dimensional 
effective potential V e ff(z) is constructed from wavefunc- 
tion of other orbitals and the equation (15) is solved for a 
new energy eigenvalue and longitudinal wavefunction. 

5.3.2. Perturbation to higher Landau levels 

With the eigenvalues and eigenfunctions of deter- 
mined, we can proceed to obtain the multi-configurational, 
single-orbital wavefunction <j> ml/ (r) using perturbation the- 
ory. We expand <p mv (r) as 

<t>mv{r) = ^ C nmv ■ Xnmvif). (18) 

n 

The coefficients c nm „ are computed by the usual non- 
degenerate perturbation theory (as is the 



„(0) 



.(1) 



„(2) 



(19) 

Matrix elements of the Oth order Hamiltonian be- 
tween different n but the same (to, v) are simply 
{nmv\h^> \n'mv) = (n 1 — n)h\OB where (r\nmv) = 
Xnmvif)- The matrix elements of become for dif- 
ferent to states due to the orthogonality of Landau func- 
tions. The perturbation terms between different v states 
are also 0, because different v states are orthogonal to each 
other. On the other hand, matrix elements of between 
different Landau states consist of two parts as follows. 

(nmu\h^ \n'mv) — (nmis\V nuc (r)\n'mis) 

+ (nmv\V eff {r)\n'mv). (20) 

For the evaluation of the matrix elements, the reader 
should refer to the appendix. 

5.4. The MCPH 3 solution for the total wavefunction and 
energy of the Hamiltonian (Type II perturbation 
theory) 

Once the Hartree single-orbital wavefunctions 4>mv(r) in 
(18) are orthonormalizcd, the Oth order total wavefunction 
is given as below since the spatial wavefunction must be 
antisymmctrized based on the FSP assumption. 



9(f) = A]l<t> mv {i f ), 



where A is the antisymmetrizing operator. 

Based on the total spatial wavefunction in (21), we 
evaluate the single-orbital energy (which is the ionization 
threshold) e m „ and total energy E to the 1st order in 
perturbation theory (called type II perturbation theory 
in §5.2). The single-orbital energy is given by, 



e - c (0) + e (i) 



(22) 



where 



\m'v' 1 1 / 

= -tE/ d3f [ d'r'^Ar^Ar)-^ 

tu'i/' 

X(t>mv{r)<Pm>v>(r')- (23) 

Only the exchange term remains in 1st order because the 
direct Coulomb term and the effective potential cancel out 
when the Hartree single-orbital wavefunctions have con- 
verged in the iteration associated with (15). 
Total energy E to 1st order is given by, 



mv 

—See 



(24) 



mv m'v' 



{^y—^rnA^v' (r ') J(25) 

In contrast, there remains part of the direct Coulomb term 
in E^ following the way we have divided the Hamiltonian 
H into H^°> and in (3). However the applicability of 

type II perturbation is determined by the ratio |eml/eml| 
as is discussed in §6.3. 



5.5. Oscillator strengths 

Oscillator strength in the length form for the transition 
AM = q is given as (M is the total magnetic quantum 
number M = J2j m -j), 



M) 
Jfi 



2m e p 2 Efi 
h~ 2 



K/l 



3 



where 



Si) 



/ 4ir 

—rjY^Oj^j). 



(26) 



(27) 



(21) 



Y[ q is the spherical harmonic. Ef t is the transition energy. 
Selection rule for the dipole-allowed transitions is : 

Am = ±1, Ai> = even 
Am — 0, Av — odd 
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5.5.1. Dipole matrices 
The dipole matrix elements in units of p are denned 



as, 



dtl, m >»> = (rn'v'\rW\mv). 



Vi r (g)i 



(28) 



\mv) and \m'v') denote an orbital in initial and final state 
respectively. Using the orthonormality of Landau func- 
tions, dipole matrix elements are evaluated as, 



(29) 



n,n' 

x / dz fm> v >{z)zfmv{z), 

+S n 'nVn + m] J dz f^ l , l/ ,(z)f mv (z), (30) 

n, n' 

+ <5„'„\/n + TO+ 1] / fm'v>( z )fmA z )- ( 31 ) 



The longitudinal wavefunction for the initial and the fi- 
nal states are not quite orthogonal to each other due to 
differing electron configurations. Nevertheless we evalu- 
ated oscillator strengths for the single electron transitions 
using only the dipole matrix element of the transitioning 
electron. The error caused by this truncation is insignifi- 
cant for one electron transitions. 

6. VALIDITY AND OTHER PHYSICAL EFFECTS 



6.1. 



Valid region in B-Z phase space 



An expansion of the single-particle Hamiltonian in terms 
of cylindrical wavefunctions provides accurate results when 
B > B c . As mentioned in §3 this is a statement that mag- 
netic field effects dominate over Coulomb field effects. On 
the other hand, the perturbative treatment of higher Lan- 
dau levels is valid when \h^\ < \h^\. Since the nuclear 
term in equation (20) dominates the effective potential 
term, the condition B > B c is also sufficient to ensure the 
validity of a perturbative treatment of the higher Landau 
levels. Therefore, B c sets a lower limit on the magnetic 
field for which our method is valid for tightly-bound elec- 
trons (i.e. in state (to, 0)). In figure 2, we plot B c for 
m = — 3 assuming Z e f f ~ Z — to for illustrative pur- 
poses. Loosely-bound electrons are well within the Landau 
regime due to their larger average distance to the nucleus. 
Therefore even at intermediate field strengths a cylindri- 
cal expansion combined with perturbative treatment of 
the higher Landau levels still provides accurate results for 
most electrons. The sharp decrease in B c with increasing 
to, seen in figure 2, means that only those tightly-bound 
inner electrons with to = — 3 will have significant config- 
uration mixing. At higher quantum state to, the Landau 
regime, with its negligible configuration mixing, is quickly 
recovered. 

6.2. Spin- flip transitions 



As the magnetic field gets smaller, spin-flip transitions 
can take place, obviating the validity of the FSP approxi- 
mation. The critical magnetic field strengths where spin- 
flip transitions take place (B s f) were recently studied for 
atoms up to Z = 10 by Ivanov & Schmclchcr (2000). In 
a fairly wide range of magnetic field they computed the 
ground state energies of atoms in which the FSP approxi- 
mation was removed - the electron spins were unrestricted. 
Using this study we found that B s f could be fit by a poly- 
nomial function of Z. At B < B s f the assumption of 
antisymmetrized total spatial wave function breaks down. 
This sets another restriction on the validity of our method 
in the phase space of (B, Z) (figure 2). 

6.3. Exchange term 

Although the significance of the exchange term in 
the binding energy increases in the FSP approximation 
(Schmelcher et al. 1999), the exchange term can be treated 
perturbatively as long as it is smaller than the 0th order 

energy eigenvalues, i.e. je^ll < \tml\- We observed this 
inequality was satisfied for any single orbital in all the 
configurations computed. In the perturbation for the ex- 
change term we did not perform calculations higher than 
0th and 1st order for the wavefunctions and energy val- 
ues because of the increasing complexity of the calcula- 
tions. Also the naive inclusion of a 1st order perturbation 
into the wavefunctions breaks the wavefunctions orthog- 
onality, requiring a further orthogonalization process for 
the single-orbital wavefunctions. The exact value of the 
wavefunction is not our primary concern, but rather the 
accuracy of the real physical observables such as energy 
values and oscillator strengths. For instance, errors on 
the energy levels caused by truncating the wavefunctions 
are actually small (§8.1). We use an explicit procedure 
to quantitatively correct the oscillator strengths in order 
to account for the deviation introduced by using Hartree 
rather than Hartree-Fock wavefunctions (§8.2). 

6.4. Electron correlation 

In reality the individual wavefunctions are not symmet- 
ric but distorted by the Coulomb coupling with other elec- 
trons. The Hartree and Hartree-Fock methods do not in- 
clude electron correlation in the Schrodingcr equation for 
each orbital. However, as the magnetic field increases, 
the electron correlation becomes less relevant (Schmelcher 
et al. 1999). The estimated error in our energies and oscil- 
lator strengths due to electron correlation is significantly 
smaller than 1% in the Landau regime. 

6.5. Relativistic effects 

The significance of relativistic effects is proportional 
to HuB/meC 2 or EB/m e c 2 where Eb is a binding en- 
ergy. The former condition translates to B/B re i (B re i = 
4.414 x 10 13 G), representing the degree of relativistic ef- 
fects in the transverse motion of electrons. However, the 
shape of the Landau wavefuntion in the relativistic theory 
is the same as in the non-relativistic theory (Lai 2000). On 
the other hand, the electron becomes relativistic along the 
magnetic field when Eb ~ m e c 2 . However, the relativistic 
effects on the Coulomb binding energies remain negligi- 
ble even at magnetic field strengths of up to 4.7 x 10 13 
G as determined by numerically solving Dirac's equation 
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for hydrogen (Lindgren & Virtamo 1979; Chen & Gold- 
man 1992). Also, the second condition is relevant only 
to the inner orbitals which arc not important for X-ray 
spectroscopy since their binding energies are high (> 10 
keV). 

6.6. Finite nuclear mass effects 

In the presence of a magnetic field, the collective motion 
of an atom and its internal degrees of freedom are coupled, 
and this coupling modifies the electronics structure of the 
atom (Pavlov & Meszaros 1993; Potekhin 1994). Conse- 
quences of this coupling include distortion of the wave- 
functions, violation of some dipole selection rules and line 
broadening. Most work on these finite nuclear mass effects 
has been performed on hydrogen atoms (Lai et al. 1992; 
Pavlov & Meszaros 1993; Potekhin 1994; Kopidakis et al. 
1996; Potekhin 1998) and hydrogenic ions (Bczchastnov 
et al. 1998). 

The immediate impact on our work is that the energy 
levels and oscillator strengths arc affected by these collec- 
tive effects. In this section we briefly review the collective 
motion and the theoretical framework for addressing it. In 
this paper we do not relax the assumption of an infinite 
nuclear mass in the Hamiltonian. This facilitates direct 
numeric comparison of our MCPH 3 approach with other 
approaches for solving the strong magnetic field Hamilto- 
nian, since they also assumed infinite nuclear mass. We 
do explain, in light of previous work on hydrogenic sys- 
tems, how we will incorporate the collective motion into 
the MCPH 3 technique in a subsequent paper. In the case 
of the bound-bound transitions of interest to us we de- 
scribe the straightforward incorporation of finite nuclear 
mass into the oscillator strength calculation. This is ap- 
propriate since previous hydrogen work also included the 
effects of the finite nuclear mass on the oscillator strengths. 

The finite nuclear mass treatment is quite complicated 
in hydrogen and has not really be dealt with previously for 
higher Z atoms. We outline in this section how substan- 
tial simplifications of this problem arise when considering 
both the specifics of higher Z systems in general and the 
requirements of X-ray spectroscopic modeling in particu- 
lar. 

6.6.1. Review of previous approaches to finite nuclear 
mass effects 

When the nuclear mass is not assumed to be infinite 
the Hamiltonian, including a term for nuclear motion, can 
be separated into center-of-mass (CM) and relative coor- 
dinates by a canonical transformation, yielding 

H —> H' = H cm + Hi n t + H c , (32) 

where H cm and Hi nt are the Hamiltonian for the CM mo- 
tion and the internal part. There is a scaling law for 
Hi n t from the infinite nuclear mass case (corresponding 
to the Hamiltonian in equation (2)) to the finite mass case 
(Pavlov- Verevkin & Zhilinskii 1980; Becken et al. 1999). 
Even when there is no translation motion of the atom, 
there is an additional term mMls (tu is the z-componcnt 
of angular momentum and £Ib is the cyclotron frequency 
of the nucleus) arising from the cyclotron motion of the nu- 
cleus (hereafter we call this term as the nuclear cyclotron 
term). Due to the nuclear cyclotron term, the large m 



states are autoionized (Kopidakis et al. 1996). Since this 
term commutes with the Hamiltonian, it is easy to imple- 
ment as an additive term in H int . Less straightforward to 
handle is the coupling term H c . This term arises because 
the motion of the CM through the magnetic field induces 
an electric field - the so-called motional Stark effect. The 
motional Stark field acts as a dipole electric field which 
increases the separation between electron and nucleus and 
therefore decreases the binding energy compared to the 
infinite nuclear mass case. H c can be expressed as : 

H c = ^- c (KxB)-f, (33) 

where M is the nuclear mass and r is the relative co- 
ordinate between the nucleus and the electron. K = 
n—e(Bxr)/c is the pseudomomentum, where n = p+eA/c 
represents the kinetic momentum. The pseudomomentum 
commutes with the Hamiltonian and therefore is a con- 
stant of motion. The pseudomomentum defines a separa- 
tion between the guiding centers of the nucleus and the 
electron; r c — -^i(K x B). The components of the pseu- 
domomentum commute for neutral atoms, while they do 
not commute for a charged system (Baye & Vincke 1990). 

The motional Stark field breaks the cylindrical symme- 
try which is present by the central Coulomb field in the 
infinite nuclear mass case. As a result the system ob- 
servables generally depend on K±, the transverse compo- 
nent of the pseudomomentum. For instance, we express 
the binding energy of an electron in an orbital (m, v) as 
e m v{K±). States can be classified according to the rela- 
tive size of compared to a critical pseudomomentum 
K c ~ (2Me ml/ (0)) 1 / 2 and e mv (0) is the binding energy of 
the electron in a non-moving atom. When K± < K c the 
central Coulomb field exceeds the Stark field; this is called 
a centered state. When K± > K c the Stark field dom- 
inates and it separates the electron from its equilibrium 
(Coulomb) position around the nucleus; this is a decen- 
tered state. 

Treatment of collective effects in the Hamiltonian de- 
pends on whether a state is centered or decentered. For 
decentered states a coordinate system f ' = r — r c is useful 
because the coupling term H c becomes zero in the trans- 
formed Hamiltonian H'. Instead the nuclear Coulomb 
term is shifted by r c , requiring evaluation of matrix ele- 
ments at each K±. For nearly centered states the motional 
Stark effect can be treated by considering the coupling 
term H c as a perturbation to the infinite nuclear mass 
solution (Pavlov & Meszaros 1993). In the perturbation 
method binding energy is given by, 

em„(#i) - e mv (0) + (34) 

Mn\ v is the transverse mass, which can be obtained by the 
2nd order perturbation, since the bound states are non- 
degenerate. The transverse mass is larger than the atomic 
mass M (mass anisotropy). 

Non-perturbative calculation was first studied by Vincke 
et al. (1992). Potekhin (1994) studied binding energies, os- 
cillator strengths and size of moving hydrogen atom by the 
multi-configurational method. Potekhin (1998) obtained 
fitting formula based on the results of Potekhin (1994). 
Bczchastnov et al. (1998) investigated hydrogen-like he- 
lium by multi-configurational method with two particle 
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basis sets. For a charged system the collective motion is 
quantized by the cyclotron motion due to the non-zero 
net charge. Therefore atomic structure is characterized by 
the Landau numbers of the collective cyclotron motion in- 
stead of K±. Baye k Vincke (1990) and Schmelcher (1995) 
have investigated a unified picture of general neutral and 
charged systems. 

6.6.2. Effect of finite nuclear mass on X-ray spectral 
features in MCPH 3 

We have only considered the infinite nuclear mass case 
in this paper to validate MCPH 3 against other methods, 
which did not consider collective effects in non-hydrogenic 
systems. However, we included a correction due to the 
nuclear cyclotron term for oscillator strengths of Am =/= 
transitions for hydrogen (mentioned elsewhere in §8.2). 
We outline here the program we will implement in a subse- 
quent paper to accommodate these collective effects within 
the MCPH 3 formalism. 

First of all, the degree of motional Stark effect decreases 
monotonically with the nuclear mass M. In addition, the 
binding energy of the inner electrons increases as the nu- 
clear charge Z becomes larger. Therefore, the tightly- 
bound states are centered or nearly centered states. This 
allows us to adopt the perturbation method for centered 
states as employed by Pavlov & Meszaros (1993) and Lai 
& Salpetcr (1995). This approach is fully compatible with 
MCPH 3 but will set a constraint on the density range of 
applicability for a given temperature. In particular this ap- 
proach is effective at higher densities where highly excited 
(decentered) states are destroyed by quasi-static electric 
fields from adjacent atoms and by collisions with fast free 
electrons (pressure ionization), or are autoionized due to 
the nuclear cyclotron term. As shown by Potekhin et al. 
(1999) the relative population of centered and decentered 
states in neutral hydrogen depends on B, p and T in a 
complicated manner. The range of validity of the pertur- 
bation theory will be the region of parameter space where 
centered states dominate. The range of parameter space in 
which a perturbative approach is valid will be much larger 
in the case of higher Z atoms compared to hydrogen be- 
cause of the heavier nuclear mass involved. We also note 
that, at sufficiently high B, differences in the nuclear cy- 
clotron term for different to states can be larger than the 
distances between single-orbital energy levels. This results 
in level mixing and sets constraint on the applicablity of 
MCPH 3 . Taking into account these points the parameter 
space (B, Z, p,T) of validity will be the subject of our next 
paper to construct realistic atmosphere models. 

A number of effects have been considered in the litera- 
ture that are much less important for X-ray spectroscopy 
in the XMM-Newton and Chandra energy bands. For in- 
stance, as the motional Stark field decreases the binding 
energy spectral features are subject to redward line broad- 
ening (magnetic broadening, Pavlov & Meszaros (1993)). 
In addition, the Stark field breaks the cylindrical sym- 
metry opening up new transitions that are forbidden ac- 
cording to the dipole selection rules for non-moving atom. 
Pavlov & Meszaros (1993) classified spectral features asso- 
ciated with the motional Stark effect as arising from high 
or low energy components. High energy components in the 
soft X-ray band consist of transitions from states of large 



binding energy; they are centered states. On the other 
hand, the low energy components of Pavlov & Meszaros 
(1993) result from decentered state transitions leading to 
radiation in the optical-UV band. Consequently we will 
be able to apply the perturbation methods to the magnetic 
broadening for spectral features in the soft X-ray band. 

7. NUMERICS 

7.1. Spatial grid points and computation time 

The main part of the iterative calculation is solving the 
eigenvalue equation (15) for each orbital and the numerical 
evaluation of the relevant integrals. We perform a finite 
element method, using a set of grid points {pi} (N p el- 
ements, grid interval Apj) and {zj} (N z elements, grid 
interval Azj) for the transverse and longitudinal direc- 
tion respectively. Since a (to, v) state is localized around 
(2m + l) 1//2 p with the width ~ p in the transverse direc- 
tion (Meszaros 1992), the transverse size of 10/5 is large 
enough for Z < 30. We precalculated all the relevant ker- 
nels (nuclear, direct and exchange) as a function of z by 
integrating over the grid points {pi} 1 since the kernels are 
independent of B when we use the cyclotron radius for 
the unit of length. See the appendix for the evaluation 
of kernels. We used a set of uniform grid points for {pi} 
with N p = 10 3 and Api = 10~ 2 . We evaluated the in- 
tegrals by a five-point Newton-Cotes integration formula. 
Numerical error associated with the kernel integrals is neg- 
ligibly small. Prior to the calculation for a given electron 
configuration, the kernels are read from the database com- 
prised of about 3x 10 5 lines of floating-point numbers. The 
kernel database is extended over m = — 30 for single- 
configurational calculation and n = 1 — 7 for the magnetic 
quantum numbers m = — 5 for multi-configurational cal- 
culation. This is sufficiently large for Z = 1 — 26 in the 
Landau regime. The algorithm of our atomic calculation 
is presented in figure 3. 

When solving an eigenvalue equation for each orbital by 
Numerov algorithm (Koonin 1986), we confined the longi- 
tudinal coordinate to < z < L, with a sufficiently large 
size of the integration box L, because f m u{z) is symmetric 
regardless of the electron correlation. Boundary condi- 
tions are : f mv (z = 0) = (Dirichlet condition) for odd 
v's, f' mv {z = 0) = for even v's (Neumann condition) and 
fmv{z = L) = for all i/'s. In contrast to the transverse 
direction, the characteristic length of longitudinal wave- 
functions varies greatly with v (Rudcrman 1971) (ao is 
the Bohr radius, a /p = a/2/3 1 / 2 ). 

L ~ (ao/ Zp) In (a /Zp) (y = 0) 

~ a v 2 /p (v > 0) 

In compensation for saving computation time from the 
kernel calculation, we confront the difficulty of the in- 
dependent size of the integration box as we keep the fine 
grid intervals. Inner tightly-bound orbitals require finer 
grid points (especially for high Z atoms), while wavefunc- 
tions of outer orbitals require large extension in integration 
box, with L as large as ~ 1000 depending on B, Z and v. 
As the results are sensitive to the grid intervals, uniform 
grid points are not appropriate for an atom mixed with 
tightly-bound and loosely-bound electrons unless we have 
a large number of grid points N gri d (meaning large com- 
putation time). Therefore, we adopt the non-uniform grid 
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points : Zj = g(zj) of Ivanov (1994). {zj} is a set of uni- 
form grid points on which the equation (15) is solved, while 
{zj} is a set of actual non- uniform grid points. A function 
g(z) must be chosen so that the grid intervals Azj become 
denser near the origin and sparser at large z. Upon the 
mapping of coordinates, the equation (15) must be trans- 
formed into an appropriate form in terms of a new set of 
grid points {zi}. For some of the loosely-bound orbitals, 
we do not know ah initio where a wavefunction changes 
steeply, but they do not require finer grid points compared 
to tightly-bound states. Therefore, we used uniform grid 
points with the interval Azj ~ 0.5 — 0.8. It reduces the 
number of grid points significantly to N z = 70 — 200 for 
tightly-bound orbitals and N z = 70 — 1000 for loosely- 
bound orbitals. Note that most electrons are usually in 
tightly-bound orbitals while a few loosely-bound orbitals 
are occupied by excited states and the ground states of 
high Z atoms. Keeping the same set of grid points for dif- 
ferent electron configurations for different B and Z is still 
a difficult task. We adjust the grid points, if necessary, 
after the 1st iteration. 

Computation time is proportional to N gr id 2 xN e 2 x Ni ter 
for the self-consistent methods due to the coupling between 
electrons (N e is the number of electrons, N it er is the num- 
ber of iterations). Use of a one-dimensional Hartree equa- 
tion is far superior to 2DHF (N grid = N p x N z = 65 x 65 = 
4225) (Ivanov & Schmclchcr 2000) in terms of the compu- 
tation time. Miller & Neuhauser (1991) adopted uniform 
grid points (N z = 500 - 1000 and Az = 0.2 - 0.8) for 
most of their electron configurations computed, suffering 
the numerical inaccuracy caused by the coarse grid points, 
as pointed out by Thurner et al. (1993). 

7.2. Convergence and numerical errors 

Our initial guess for the longitudinal wavefunction of 
the (to, v) orbital is adopted from restricted variational 
studies (Flowers et al. 1977). 

f mv {z) ex z"e- a -^ (35) 

The final results do not change with our initial choice of 
{a>mv} varying from 0.1 to 10. Energy values converge to 
less than 1% after the 2nd iteration and less than 0.1% 
after the 4th iteration, although convergence speed be- 
comes slow as N e increases (figure 4). Therefore, estima- 
tions of convergence speed from hydrogen results (Miller 
& Neuhauser 1991) arc not correct. Convergence speed 
for a 2DHF calculation will be slower than in the one- 
dimensional Hartree method. 

Upon the evaluation of the direct and exchange terms, 
we found that the coupling terms between different n > 
states are negligible (< 0.1%). Therefore, we did not in- 
clude the contribution from n > Landau states of other 
orbitals (setting n = n' = in (A9) and (A10) in the ap- 
pendix), reducing the computation time when evaluating 
the 1st order total energy. The numerical error is overall 
less than 0.1%. 

8. RESULTS 

8.1. Energy values and ground state configurations 

We have computed total energies for ground states and 
several excited states at representative magnetic fields. 
The computed energies match with the previously most 



accurate values (e.g. 2DHF) with less than 1% deviation. 
In general, model-dependent total energy is closer to the 
exact value when it is lower than other models. However, 
the computed energy values can be lower than the exact 
values due to numerical error or inconsistency in atomic 
models. At pz close to 1, our energy values are actually 
lower than the previous results known to be accurate to 4 
digits. This is because higher Landau states of the inner- 
most electrons are no longer perturbations on the n = 
state. Therefore, perturbative method may provide total 
energies lower than the previous accurate data. Neverthe- 
less, we listed the multi-configurational results as long as 
the discrepancy is less than 1%. 

8.1.1. Hydrogen 

We present the energies for ground states and several ex- 
cited states in comparison with the multi-configurational 
Hartree-Fock method (MCHF hereafter) (Ruder et al. 
1994) (table 1, 2). They fully included higher Landau 
levels up to n = 12 in the Hartree-Fock equation. We pro- 
duced fairly good energy values with < 1% deviation from 
MCHF. 

8.1.2. Helium 

We compare our results with 1DHF (Thurner et al. 
1993; Miller & Neuhauser 1991) and UHF (Jones et al. 
1999b) for the ground state and several excited states (ta- 
ble 3, 4). The ground state configuration for helium is 
{m,v) — (0,0), (1,0) in the Landau regime. Our single- 
configurational total energies (excluding n > Landau 
levels) match well with 1DHF results (Thurner et al. 1993) 
by < 0.5%. Our multi-configurational total energies are 
lower than 1DHF by <~ 1 — 5% at f3z < 10 and higher than 
UHF by < 1% at Z > 1. 

8.1.3. High Z atoms 

For Z > 2 atoms, published atomic data is scarce 
(Neuhauser et al. 1987; Ivanov & Schmelcher 2000). First 
of all, we searched for the ground state configurations of 
neutral atoms Z = 2 — 26 at some representative field 
strengths [table 5,6]. In the limit of B — > oo, the ground 
state configuration is such that all the bound electrons are 
in tightly-bound orbitals (m,v) = (0, 0), (1, 0)...(N e - 1, 0). 
With decreasing B, the spatial crossover takes place, i.e. 
having electrons in loosely-bound orbitals rather than in 
outer tightly-bound orbitals, lowers the total energy. A 
v = 2 orbital may form the ground state configuration 
because a v — 2 orbital can be closer to the nucleus than 
v = 1 orbitals due to the electron screening (Jones 1985b). 
However, we did not find any ground state configurations 
with v = 2 orbitals at the magnetic fields in tables 5 and 
6. Several spatial crossovers take place before the spin-flip 
transition sets in, as we observed for high Z atoms (Ivanov 
& Schmelcher 2000). Ivanov & Schmelcher (2000) investi- 
gated the ground state configurations at (3 <~ 10 _1 — 10 4 
for Z = 1 — 10 atoms. For Z > 10 atoms, only Neuhauser 
et al. (1987) provided the ground state configurations and 
the total energies among the self-consistent methods. 

Our single-configurational total energies match very well 
with 1DHF (Neuhauser et al. 1987) up to Z = 5, while 
1DHF total energies become higher than our values by 
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at most 1% at Z > 5. This is probably because their grid 
points are not sufficiently fine for inner electrons in high Z 
atoms. Our multi-configurational results are 5 — 10% lower 
than 1DHF at (3 Z < 10 and a few % lower at (3 Z < 100. 
Some of their ground state energies are even higher than 
our multi-configurational energy of several excited states. 

Jones (1986) studied the ground state energy for several 
high Z neutral atoms by DF. Their results are a few % 
lower than our results and 2DHF for the ground states of 
several atoms, while they fixed n = (table 6). Ground 
state energies obtained by DF are highly dependent on 
the form of the exchange-correlation density functionals 
(Jones 1985a), even though DF can relatively easily im- 
plement the exchange and electron correlation terms in 
the local density approximation. 

Up to Z — 10, our multi-configurational ground state 
energies are higher than 2DHF results by < 1% overall 
at (3 'z > I- Even though 2DHF claims no numerical er- 
rors by the use of the Richardson extrapolation method, 
this difference could be attributed to the fact that we did 
not include the 2nd order perturbation in the total en- 
ergy. The errors are expected to come mainly from the 
inner electrons but to be smaller for outer electrons. How- 
ever, the inner electrons of high Z atoms are irrelevant to 
X-ray spectroscopy since their binding energy is relatively 
high (> 10 keV), while the inner electrons in low Z atoms 
can contribute to the opacities in the X-ray band. The 
accuracy in total energies within 1% in comparison with 
2DHF results is sufficient for the identification of the po- 
tential atomic lines or edges. Our energies as well as the 
0th order wavefunctions can be improved by adding an ap- 
propriate local exchange term to a single-orbital Hartree 
equation. In contrast, empirical density functionals for 
the exchange term in the Landau regime are not known 
presently. 

8.2. Oscillator strengths 

Compared to energy values, oscillator strengths are sub- 
ject to more atomic model dependent uncertainties. This 
is because atomic models are usually optimized for mini- 
mizing the total energy The completeness of the orbital 
wavefunctions represented by Thomus-Kuhn rules is a sec- 
ondary issue. Errors in transition energy values can affect 
the oscillator strengths. For the candidate strong transi- 
tions, 10% accuracy in the oscillator strengths is sufficient, 
considering that spectral features can undergo a number 
of broadening effects (Rajagopal et al. 1997). For the same 
reasons, we do not consider transitions other than dipole- 
allowed. For hydrogen, we compare to the results from 
MCHF (Ruder ct al. 1994) and we included a correction 
by the proton cyclotron term for Am ^ transitions. Con- 
sequently, the oscillator strengths match fairly well with 
their results (table 7). 

Data for oscillator strengths at Z > 1 were found only 
for 1DHF. Miller (1990) published data at Z = 2, 6 for 
bound-bound transitions, while Ruder et al. (1994) ob- 
tained the oscillator strengths only for helium-like atoms. 
It turns out that tight-tight transitions match with 1DHF 
very well, while the oscillator strengths for tight-loose 
transitions are consistently larger than 1DHF results by 
10 — 50 %. The most likely explanation is that we did 



not include the exchange effects in the orbital wavefunc- 
tion. In the antisymmetrized Hartree total wavefunctions, 
electrons are closer to each other than in the Hartree-Fock 
total wavefunctions. As exchange effects are larger for a 
pair of closer electrons, tightly-bound electrons are more 
deviant from Hartree-Fock orbitals. This discrepancy is 
small for the inner tightly-bound electrons since their be- 
havior is mainly determined by the Coulomb attraction 
from the nucleus. Also, it is very small for a loosely-bound 
electron simply because it is far from other electrons, un- 
less there are other loosely-bound electrons nearby. For 
tight-tight transitions, a transitioning electron (usually the 
outermost tightly-bound electron) in the Hartree scheme 
is closer to the nucleus than in the Hartree-Fock scheme 
both in the initial and final state. This deviation can- 
cels out when we evaluate the oscillator strengths by the 
length form because a change in the spatial distribution 
of the electron associated with the atomic transition is 
about the same in the Hartree and the Hartree-Fock pic- 
ture. However, for tight-loose transitions, the deviation in 
the initial state remains in the oscillator strengths and it 
gives us the consistently larger oscillator strengths because 
the Hartree scheme requires larger change in the spatial 
distribution compared to the Hartree-Fock scheme. 

A well-defined measure of this deviation associated with 
each bound electron is the ratio of the exchange term to 
the Hartree single-orbital energy, i.e. r ex = \ern U /^mL\- 
Therefore, we introduced a correction to the oscillator 
strengths for the tight-loose transitions by the amount of 
r ex . The same correction is applicable when size of atoms 
or ions is required to include non-ideal effects, for instance 
through the excluded- volume term (Potckhin et al. 1999). 
In contrast to hydrogen case, we did not include the fi- 
nite nuclear mass correction since the previous results for 
comparison assumed the infinite nuclear mass. As a re- 
sult, the corrected oscillator strengths match with 1DHF 
results within 10% for different tight-loose transitions of 
helium and carbon at various magnetic fields (table 8, 9). 
The discrepancy at small 0z is due to the finite occu- 
pation in higher Landau levels, which is not included in 
1DHF but is included in MCPH 3 . Although comparisons 
are available only for helium and carbon due to few pub- 
lications on oscillator strengths for Z > 1, the numerical 
results in tables 8 and 9 based on this physically reason- 
able correction, strongly support that oscillator strengths 
are calculated with 10% accuracy for other atoms up to 
Fe. 

9. FUTURE WORK 

Construction of an extensive atomic data base from the 
method presented in this article is in progress. It consists 
of energies and oscillator strengths for a large number of 
grid points in B. Based on the atomic data, equation of 
state and opacity tables will be constructed in the phase 
space of (B, p, T, Z) covering the typical parameter values 
for cooling neutron stars. 

We are grateful to Ehud Behar for careful reading the 
manuscript. This work was partially supported by NASA 
grant NAG5-7737. 
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APPENDIX 

EVALUATION OF POTENTIALS AND ENERGIES 

Analytical expression of potentials and energies 
In order to avoid the complications caused by the r _1 cusp at the origin, we expanded the terms r _1 and \ f— r' l_l 



appearing in the potentials by the Bessel functions J m (x) (Flowers et al. 1977). L m (x), L™ (i) denote the Laguerre function 
and associated Laguerre function respectively. Hereafter, J dk = J °° dk, J dz = dz, J dp — J °° dp, J d<j> — J Q V dcj> . 

L= = fdkJ (kp)e- k \'\, (Al) 
+ z z J 



r- 1 



Vp 2 



If-r'r 1 



yjp 2 + p' 2 - 2pp' cos(<P - <//) + (z- z'f 

OO „ 

I dke m ^'^J m {kp)J m {kp')e- k \ z - z '\. (A2) 



m— — oo 



For the single-configurational potentials, 



Ze 2 

V nuc (z) = yV mQ0 {z), (A3) 

V °ffW = iE I dz ' D ^oo(z ~ z')\f m , v ,{z')\ 2 , (A4) 



where nuclear and direct kernels are defined as 



V mnn ,{z) = J pdp J d4^* n , m {p,<j>)^ nm {p,<j>) (A5) 
D^ w (z-z')= J pdp Jd<t>J p'dp'J d<j>' $>l ml {p',<p')3> Wm ,{p',<p') 

X$» m (p, <f>)K>m(P, ^ l^^^l ' (A6) 

The matrix elements for the multi-configurational calculation are : 

Ze 2 f 

{nmv\V nuc (r)\n'mv) = — / dz V mnn >{z)\f mv {z)\ 2 , (A7) 

P J 

{nmv\V efs {r)\rimv) = % ^ ^ Chm'v'Cn'm'v' j dz D™?£^,{z - z') 

X\fm v (z)\ 2 \fm^(z')\ 2 . (AS) 

1st order perturbation of single-orbital energy and total energy are given as, 



^rnv - ^ ^ } ^ ^ ^ Cnmi/Cn'mi/Cnrn'i/'Cn'rn'v'' J d,Z J dz E^^^/^Z Z ) 

^ rri'v' n,n' h,h' 

Xfm V {z)f m 'u'{z')f^, v ,(z)f^ v (z'), (A9) 

- _£l W W c c / c- , ,c-> > > I dz I dz' 

mu m'v' n,n' h,h' 



yD^Uz~ z')\f mv {z)\ 2 \f m , v> {z')\ 2 

+E™ n i,(z - z')f mu {z)f mlv ,(z')r m ,Az)r m Az')^ , (A10) 
where exchange kernels are defined as 

EZ^{z-z') ee J pdp I d4> I p'dp' I d<j>' *', m (J,4/)*i, m ,(f/,<t/) 

X p l — s $nm{p,4>)$n m '{p,4>)- (All) 

r — r ' 



12 




13 



REFERENCES 



Abrahams, A. M. & Shapiro, S. L. 1991, ApJ, 382, 233 
Alcock, C. & Illarionov, A. 1980, ApJ, 235, 534 
Baye, D. & Vincke, M. 1990, Phys. Rev. A, 42, 391 
Beckcn, W. & Schmelcher, P. 2000, J. Phys. B, 33, 545 
— . 2001, Phys. Rev. A, 63, 0a4+ 

Becken, W., Schmelcher, P., & Diakonos, F. K. 1999, J. Phys. B, 
32, 1557 

Becker, W. & Triimper, J. 1997, A&A, 326, 682 
Bczchastnov, V. G., Pavlov, G. G., & Ventura, J. 1998, 

Phys. Rev. A, 58, 180 
Camilo, F., Thorsett, S. E., & Kulkarni, S. R. 1994, ApJ, 421, L15 
Chen, H., Ruderman, M. A., & Sutherland, P. G. 1974, ApJ, 191, 

473 

Chen, Z. & Goldman, S. P. 1992, Phys. Rev. A, 45, 1722 
Chevalier, R. A. 1996, ApJ, 459, 322 

Demeur, M., Heenen, P. ., & Godefroid, M. 1994, Phys. Rev. A, 49, 
176 

Flowers, E. G., Ruderman, M. A., Lee, J. ., Sutherland, P. G., 
Hillebrandt, W., & Mueller, E. 1977, ApJ, 215, 291 

Forster, H., Strupat, W., Rosner, W., Wunner, G., Ruder, H., & 
Herold, H. 1984, J. Phys. B, 17, 1301 

Glasscr, M. L. 1975, ApJ, 199, 206 

Glasscr, M. L. & Kaplan, J. I. 1975, ApJ, 199, 208 

Hewish, A. S., Bell, A. J., Pilkington, J. D. H., Scott, P. F., & 
Collins, R. A. 1968, Nature, 217, 709 

Ho, W. C. G. & Lai, D. 2001, preprint (astro-ph/0104199) 

Holas, A. & March, N. H. 1997, Phys. Rev. A, 56, 4595 

Hylleraas, E. A. 1930, Z. Physik, 65, 209 

Ivanov, M. V. 1994, J. Phys. B, 27, 4513 

Ivanov, M. V. & Schmelcher, P. 1998, Phys. Rev. A, 57, 3793 
— . 1999, Phys. Rev. A, 60, 3558 
— . 2000, Phys. Rev. A, 61, 022505 

Johnscn, K. & Yngvason, J. 1996, Phys. Rev. A, 54, 1936 
Jones, M. D., Ortiz, G., & Ceperley, D. M. 1997, Phys. Rev. E, 55, 
6202 

— . 1999a, A&A, 343, L91 

— . 1999b, Phys. Rev. A, 59, 2875 

Jones, P. B. 1985a, Phys. Rev. Lett., 55, 1338 

— . 1985b, MNRAS, 216, 503 

— . 1986, MNRAS, 218, 477 

Jordan, S., Schmelcher, P., Becken, W., & Schweizer, W. 1998, 
A&A, 336, L33 

Koonin, S. E. 1986, Computational Physics (Bcnjamin/Cummings, 

Menlo Park, California) 
Kopidakis, N., Ventura, J., & Herold, H. 1996, A&A, 308, 747 
Kossl, D., Wolff, R. G., Miiller, E., & Hillebrandt, W. 1988, A&A, 

205, 347 

Lai, D. 2000, preprint (astro-ph/0009333) 

Lai, D. & Salpetcr, E. E. 1995, Phys. Rev. A, 52, 2611 

— . 1996, Phys. Rev. A, 53, 152 

— . 1997, ApJ, 491, 270 

Lai, D., Salpeter, E. E., & Shapiro, S. L. 1992, Phys. Rev. A, 45, 
4832 

Landau, L. D. & Lifshitz, E. M. 1965, Quantum mechanics (Course 
of theoretical physics, Oxford: Pergamon Press, 1965) 

Lieb, E. H., Solovej, J. P., & Yngvason, J. 1992, Phys. Rev. Lett., 
69 749 

Lindgren, K. A. U. & Virtamo, J. T. 1979, J. Phys. B, 12, 3465 
Meszaros, P. 1992, High-energy radiation from magnetized neutron 

stars (Theoretical Astrophysics, Chicago: University of Chicago 

Press, — cl992) 

Miller, M. C. 1990, PhD thesis, California Inst, of Tech., Pasadena. 
— . 1992, MNRAS, 255, 129 

Miller, M. C. & Ncuhauscr, D. 1991, MNRAS, 253, 107 
Miiller, E. 1984, A&A, 130, 415 

Ncuhauscr, D., Koonin, S. E., & Langankc, K. 1987, Phys. Rev. A, 
36, 4163 



Ncuhauscr, D., Langanke, K., & Koonin, S. E. 1986, Phys. Rev. A, 
33, 2084 

Ortiz, G., Jones, M. D., & Ceperley, D. M. 1995, Phys. Rev. A, 52, 
3405 

Ostriker, J. P. & Gunn, J. E. 1969, Nature, 221, 454 
Ozel, F. 2001, preprint (astro-ph/0103227) 

Paerels, F., Mori, K., Motch, C, Haberl, F., Zavlin, V. E., Zane, S., 
Ramsay, G., Cropper, M., & Brinkman, B. 2001, A&A, 365, L298 

Page, D. 1998, in NATO ASIC Proc. 515: The Many Faces of 
Neutron Stars., 539 

Pavlov, G. G. & Meszaros, P. 1993, ApJ, 416, 752 

Pavlov, G. G., Shibanov, Y. A., Zavlin, V. E., & Meyer, R. D. 
1995, in The Lives of the Neutron Stars. Proceedings of the 
NATO Advanced Study Institute on the Lives of the Neutron 
Stars, Publisher, Kluwer Academic, Dordrecht, The Netherlands, 
Boston, Massachusetts, 1995. p. 71 

Pavlov, G. G., Zavlin, V. E., Sanwal, D., Burwitz, V., & Garmirc, 
G. 2001, ApJ, 552, L129 

Pavlov- Vercvkin, V. B. & Zhilinskii, B. 1980, Phys. Lett. A, 78, 244 

Pons, J. A., Walter, F. W., Lattimcr, J. M., & Prakash, M. 2000, 
AAS/High Energy Astrophysics Division, 32, 3314 

Potekhin, A. Y. 1994, J. Phys. B, 27, 1073 

— . 1998, J. Phys. B, 31, 49 

Potekhin, A. Y., Chabrier, G., & Shibanov, Y. A. 1999, 

Phys. Rev. E, 60, 2193 
Rajagopal, M. & Romani, R. W. 1996, ApJ, 461, 327 
Rajagopal, M., Romani, R. W., & Miller, M. C. 1997, ApJ, 479, 347 
Rclovsky, B. M. & Ruder, H. 1996, Phys. Rev. A, 53, 4068 
Rognvaldsson, O. E., Fushiki, I., Gudmundsson, E. H., Pethick, 

C. J., & Yngvason, J. 1993, ApJ, 416, 276 
Romani, R. W. 1987, ApJ, 313, 718 

Rosner, W., Wunner, G., Herold, H., & Ruder, H. 1984, J. Phys. B, 
17, 29 

Ruder, H., Wunner, G., Herold, H., & Geyer, F. 1994, Atoms in 
Strong Magnetic Fields. Quantum Mechanical Treatment and 
Applications in Astrophysics and Quantum Chaos (X, 309 pp. 93 
figs.. Springer- Vcrlag Berlin Heidelberg New York.) 

Ruderman, M. A. 1971, Phys. Rev. Lett., 27, 1306 

Salpeter, E. E. 1998, Journal dc Physique, 10, 11285 

Schmelcher, P. 1995, Phys. Rev. A, 52, 130 

Schmelcher, P., Ivanov, M. V., & Becken, W. 1999, Phys. Rev. A, 
59, 3424 

Scrinzi, A. 1998, Phys. Rev. A, 58, 3879 
Skjcrvold, J. E. & Ostgaard, E. 1984, Phys. Scr, 29, 543 
Steinberg, M., Ortner, J., & Ebeling, W. 1998, Phys. Rev. E, 58, 
3806 

Thorolfsson, A., Rognvaldsson, O. E., Yngvason, J., & 

Gudmundsson, E. H. 1998, ApJ, 502, 847 
Thurner, G., Korbel, H., Braun, M., Herold, H., Ruder, H., & 

Wunner, G. 1993, J. Phys. B, 26, 4719 
Triimper, J., Pietsch, W., Rcppin, C, Voges, W., Staubert, R., & 

Kcndziorra, E. 1978, ApJ, 219, L105 
Tsuruta, S. 1995, in The Lives of the Neutron Stars. Proceedings of 

the NATO Advanced Study Institute on the Lives of the Neutron 

Stars, Publisher, Kluwer Academic, Dordrecht, The Netherlands, 

Boston, Massachusetts, 1995. p. 133 
Vincke, M., LeDourneuf, M., & Baye, D. 1992, Journal of Physics B 

Atomic Molecular Physics, 25, 2787 
Wheaton, W. A., Doty, J. P., Primini, F. A., Cooke, B. A., Dobson, 

C. A., Goldman, A., Hecht, M., Howe, S. K., Hoffman, J. A., & 

Scheepmaker, A. 1979, Nature, 282, 240 
Woosley, S. E. & Weaver, T. A. 1986, ARA&A, 24, 205 
Yakovlcv, D. G., Levenfish, K. P., & Shibanov, Y. A. 1999, 

Physics-Uspekhi, 42, 737 
Zane, S., Turolla, R., & Treves, A. 2000, ApJ, 537, 387 
Zavlin, V. E. & Pavlov, G. G. 1998, A&A, 329, 583 
Zavlin, V. E., Pavlov, G. G., & Shibanov, Y. A. 1996, A&A, 315, 

141 



14 



Fig. 1. — Comparison of various atomic models in the phase diagram of Z and Pz- 1DHF = Hartree-Fock method with the adiabatic 
approximation (Neuhauser et al. 1986). 2DHF = Two-dimensional Hartree-Fock method (Ivanov &; Schmelcher 2000). DF = Density 
functional method. TFD = Thomas-Fermi-Dirac method. 
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Fig. 2. — Critical magnetic field B c for a tightly-bound orbital (m, 0) and the spin-flip transition field strength B a f. We assumed 
Z e ff ~ Z — m for illustrative purpose. We did a polynomial fit to the spin-flip field strengths obtained by Ivanov & Schmelcher (2000). 
B sf = 1.7 X 10 9 G for Z = 2, B sf = 9.6 X 10 10 G for Z = 10 and B s} = 6.8 X 10 11 G for Z = 26. Triangle points represent B sf for neutral 
atoms, while crossed points are for singly-ionized atoms. 
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Table 1 

Energies of tightly-bound states for hydrogen. 



(m, v) 


(0,0) 




(1,0) 




(2,0) 




(3,0) 




/3 


Present work 


R.udcr a 


Present work 


Ruder 


Present work 


Ruder 


Present work 


Ruder 


1 


28.10 


27.82 b 


16.33 


16.32 b 


12.81 


12.82 b 


10.84 


10.97 b 


2 


34.53 


34.85 b 


21.38 


21.44 b 


16.84 


17.06 b 


14.57 


14.69 


5 


47.00 


47.56 b 


30.49 


30.62 b 


24.52 


24.72 


21.33 


21.44 


10 


59.61 


60.26 b 


39.49 


39.88 


32.31 


32.48 


28.22 


28.32 


20 


75.50 


76.22 b 


51.26 


51.60 


42.24 


42.39 


37.06 


37.14 


50 


102.41 


103.11 


71.38 


71.70 


59.42 


59.54 


52.44 


52.52 


100 


127.93 


128.64 


90.74 


91.08 


76.13 


76.25 


67.52 


67.59 


200 


158.50 


159.24 


114.45 


114.70 


96.66 


96.77 


86.14 


86.20 


500 


207.47 


208.52 


153.17 


153.43 


130.62 


130.76 


117.15 


117.24 


1000 


252.30 


253.21 


188.72 


189.18 


162.17 


162.42 


146.17 


146.34 



a Multi-configurational Hartrcc-Fock method (Ruder ct al. 1994) 
b Encrgies are calculated by the spherical expansion of wave functions. 



Note. — All energies are in cV. A.z — 2 X 10 



15 



FlG. 3. — Algorithm of the Multi-configurational Perturbative Hybrid Hartree Hartrec-Fock method. 



Initial setup 

1. Grid points 

2. Trial function 



A given set of atoms, 
electron configurations 
and magnetic field 




Kernel storage 

1. Nuclear kernels 

2. Direct kernels 

3. Exchange kernels 



1. Compute effective potentials 

2. Solve eigenvalue equations 




Type I perturbation 



Configuration mixing 
with higher Landau levels 



Antisymmetrization 



Total wavefunction 



Total energy 



Type II perturbation 
"I 



Output data 

1. Transition & ionization energies 

2. Oscillator strengths 



FlG. 4. — Convergence test for the single-orbital Hartree energy of the innermost orbital (m, v) = (0, 0) for Z = 2, 6, 14, 26 neutral atoms. 
B12 = 1, Az = 4 X 10 -2 . All the electrons are in tightly-bound orbitals. There was minor change in the convergence speed for different 
magnetic fields and electron configurations. 



6x° 
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Table 2 

Energies of loosely-bound states for hydrogen. 



(m.u) (0,1) (0,2) (1,1) (1,2) 






Present work 


Ruder a 


Present work 


Ruder 


Present work 


Ruder 


Present work 


Ruder 


1 


8.051 


8.101 b 


4.584 


4.733 b 


6.653 


6.673 


3.832 


3.876 


2 


9.096 


9.135 b 


5.053 


5.138 b 


7.759 


7.778 


4.314 


4.346 


5 


10.38 


10.41 


5.650 


5.684 


9.221 


9.224 


4.940 


4.961 


10 


11.24 


11.24 


6.076 


6.090 


10.23 


10.23 


5.396 


5.411 


20 


11.94 


11.94 


6.480 


6.481 


11.10 


11.12 


5.834 


5.845 


50 


12.62 


12.62 


6.956 


6.971 


12.05 


12.05 


6.381 


6.388 


100 


12.96 


12.96 


7.310 


7.319 


12.55 


12.56 


6.767 


6.775 


200 


13.20 


13.20 


7.632 


7.647 


12.94 


12.94 


7.136 


7.140 


500 


13.39 


13.40 


8.040 


8.051 


13.24 


13.25 


7.590 


7.590 


1000 


13.49 


13.49 


8.329 


8.333 


13.39 


13.39 


7.907 


7.906 



a Multi-configurational Hartrcc-Fock method (Ruder ct al. 1994) 
b Encrgics arc calculated by the spherical expansion of wave functions. 



Note. — All energies are in eV. Az — 0.2. 



Table 3 



Ground state energies for helium. 



0z 


Present work 


Ruder a 


Jones b 


1 


146.06 (134.36) 


134.29 


146.94 


2 


184.73 (176.34) 


176.30 


186.85 


5 


255.39 (249.71) 


249.80 


258.19 


10 


326.05 (322.31) 


322.22 


329.27 


20 


414.92 (412.19) 


412.26 


418.22 


50 


564.50 (562.74) 


563.61 


568.39 


100 


707.15 (705.91) 


706.99 


710.95 


200 


880.00 (879.14) 


879.04 




500 


1158.49 (1158.49) 


1156.63 




1000 


1410.90 (1410.90) 


1409.13 





a lDHF (Ruder ct al. 1994) 
b UHF (Jones et al. 1999b) 



Note. — All energies are in eV. The ground 
state configuration is (m, v) = (0,0), (1,0). 
Pz = 1 f° r helium corresponds to 1.88 x 10 10 
G. Bracketed values are results from our singlc- 
configurational calculation. Az = 4 x 10~ 2 . 

Niter = 3. 



Table 4 



Energies of excited states for helium at j3z = 50. 



(m,z/) 


Present work 


Miller a 


Ruder b 


Jones c 


(0,1) 


423.29 


413.7 


421.71 


426.39 


(0,2) 


418.34 


408.8 


416.17 


420.88 


(1,1) 


421.72 


412.2 


419.68 


424.38 


(1,2) 


416.61 


407.8 


415.25 


419.95 


(2,0) 


521.73 




521.37 


525.97 


(2,1) 


421.45 




419.49 


424.18 


(2,2) 


416.26 




414.83 


419.52 



a lDHF (Miller & Neuhauser 1991) 
b lDHF (Ruder et al. 1994) 
C UHF (Jones ct al. 1999b) 



Note. — All energies are in eV. One electron is in 
(m,z/) = (0,0). AS = 4 x 10~ 2 . N iter = 3. 



Table 5 

Ground state energies of neutral atoms through Z = 1 — 14. 







B12 = 0.1 






Bl2 


= 0.5 




z 


Present work 


Ivanov a 


Neuhauser b 


Present work 


Ivanov 


Neuhauser 


1 


0.0770 (0.0762) 


0.07781 


0.0761 


0.1301 


(0.1298) 


0.13114 


0.130 


2 


0.2613 (0.2557) 


0.26387 


0.255 


0.4567 


(0.4543) 


0.46063 


0.454 


3 


0.5337 (0.5160) 


0.54042 


0.516 


0.9526 


(0.9449) 


0.96180 


0.944 


4 


0.8897 (0.8468) 


0.89833 


0.846 


1.600 


(1.580) 


1.61624 


1.580 


5 


1.327 (1.238) 


1.33229 


1.238 


2.390 


(2.349) 


2.41101 


2.347 


6 


1.854 (1.687) 


1.83895 


1.678 


3.308 


(3.237) 


3.33639 


3.22 


7 


(2.186) 


2.41607 


2.17 


4.353 


(4.239) 


4.38483 


4.22 


8 


(2.754 [1]) 


3.08253 [1] 


2.71 


5.517 


(5.338) 


5.55032 


5.32 


9 


(3.377 [1]) 


3.82966 [1] 


3.36 


6.803 


(6.544) 


6.82794 


6.51 


10 


(4.049 [1]) 


4.65087 [1] 




8.198 


(7.834) 


8.21365 


7.819 


11 








9.718 


(9.224) 




9.197 


12 








11.410 [1] 


(10.750 [1]) 




10.72 [1] 


13 








13.251 [1] 


(12.364 [1]) 




12.32 [1] 


14 








15.246 [2] 


(14.062 [2]) 




14.00 [1] 



a 2DHF (Ivanov & Schmelcher 2000) 
b lDHF (Neuhauser et al. 1986) 



Note. — All energies are in keV. The values in parentheses are single-configurational energy values. 
The number of v = 1 orbitals are in square brackets (No brackets indicate all the electrons are in 
v = Q orbital). A5 = (2 - 5) x 10~ 2 . N it er = 5. 
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Table 6 

Ground state energies of neutral atoms through Z = 1 — 26. 









B±2 — 


1.0 










Bl2 = 


5.0 




z 


Present work Ivanov a 


Ncuhauscr b Jones c 


Present work Ivanov 


Neuhauscr 


Jones 


1 


0.1614 


0.16222 


0.161 




0.2558 


0.25750 


0.255 




2 


0.5766 


0.57999 


0.574 




0.9574 


0.96191 


0.958 


1.040 


3 


1.214 




1.22443 


1.209 




2.078 




2.08931 


2.076 




4 


2.056 




2.07309 


2.042 




3.586 




3.61033 


3.584 




5 


3.085 




3.10924 


3.054 




5.476 




5.49950 


5.456 




6 


4.288 




4.31991 


4.20 






7.695 




7.73528 


7.60 


8.03 


7 


5.657 




5.69465 


5.54 






10.231 


10.29919 


10.20 




8 


7.176 




7.22492 


7.02 






13.099 


13.17543 


13.00 




9 


8.845 




8.90360 


8.63 






16.264 


16.34997 


16.10 




10 


10.664 


10.72452 


10.39 


10.70 


19.702 


19.81072 


19.57 


20.24 


11 


12.625 




12.25 




23.406 




24.64 




12 


14.745 




14.23 




27.436 




27.17 




13 


16.973 


1 




16.34 


1 




31.675 




31.35 




14 


19.408 


1 




18.60 


1 


19.09 


36.154 




35.74 


36.76 


15 


21.987 


1 




20.95 


1 




40.915 




40.35 




16 


24.718 


2 




23.43 


2 




45.881 




45.22 




17 


27.618 


2 




26.07 


2 




51.067 




50.30 




18 


30.766 


2 




28.82 


2 




56.530 




55.95 




19 


34.036 


2 










62.181 








20 


37.500 


? 








35.48 


68.031 






68.37 


21 












74.184 


1 








22 














80.602 


1 








23 














87.263 


1 








24 














94.259 


1 








25 














101.25 


1 


... 






26 


(55.410 [5]) 


55.10 [6 


56.01 


108.64 


2 


... 


106.09 [2] 


108.18 



a 2DHF (Ivanov & Schmelcher 2000) 
b lDHF (Neuhauser et al. 1986) 
°DF (Jones 1985b) 



Note. — All energies arc in kcV. Az — (2 — 5) X 10 2 . Nu er — 5. Other atomic models computed several 
ground state energies. (1) Thomas- Fcrmi-Dirac method (Skjcrvold & Ostgaard 1984): 4.14, 7.73 [kcV] for Z — 6 
at B12 - 1.0,5.0. 56.21, 105.89 [keV] for Z - 26 at B 12 = 1.0,5.0. (2) Restricted variational method (Flowers 
ct al. 1977): 0.545, 0.913 [keV] for Z = 2 at B 12 = 1.0,5.0. 53.13, 101.7 [keV] for Z = 26 at B 12 = 1.0,5.0. (3) 
Density functional method with correction (Kossl ct al. 1988) provide lower energy value 110.4 [keV] for Z — 26 
at B12 = 5.0. 



Table 7 

Oscillator strengths of bound-bound transitions for hydrogen. 



(m, v) — > (m' , v') 


P = 5 


P = 10 


(3 = 50 


(3 = 100 


P = 500 


(0,0) 


- (1,0) 


0.115 


7.15xl0~ 2 


2.29xl0" 2 


1.46xl0~ 2 


5.04xl0~ 3 






(0.108) 


(6.88xl0~ 2 ) 


(2.32xl0~ 2 ) 


(1.44xl0~ 2 ) 


(5.06xl0~ 3 ) 


(0,0) 


- (0,1) 


0.627 


0.553 


0.437 


0.364 


0.228 


(0.615) 


(0.560) 


(0.414) 


(0.351) 


(0.227) 


(0,0) 


-> (0,3) 


5.28xl0~ 2 


4.93xl0~ 2 


4.19xl0~ 2 


3.64xl0~ 2 


2.47xl0~ 3 






(5.24xl0~ 2 ) 


(4.94xl0~ 2 ) 


(3.99xl0~ 2 ) 


(3.52xl0~ 2 ) 


(2.45xl0~ 3 ) 


(1,0) 


- (2,0) 


8.80xl0~ 2 


5.47xl0~ 2 


1.97xl0~ 2 


1.30xl0~ 2 


5.47xl0~ 3 






(8.46xl0~ 2 ) 


(5.47xl0~ 2 ) 


(1.98xl0~ 2 ) 


(1.30xl0~ 2 ) 


(5.48xl0~ 3 ) 


(0,1) 


- (0,2) 


1.53 


1.50 


1.38 


1.33 


1.20 






(1.52) 


(1.48) 


(1.38) 


(1.33) 


(1.20) 



Note. — Grid points are the same as in the tables of energy values. The values in parentheses 
are from multi-configurational Hartree-Fock method (Ruder et al. 1994). For Am = 1 transitions, 
finite nuclear mass correction is included. 
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Table 8 

Oscillator strengths of bound-bound transitions for helium. 



(to, 


") 


/?z = 5 


Pz = 10 


Pz = 50 


pz = 100 


Pz = 500 


(1.0) - 


- (1,1) 


0.523 


0.424 


0.240 


0.187 


0.103 


CD 5001 


(0.410) 


CO 2391 


(0.184) 


CO 0991 


(1.0) - 


- (2,0) 


7.52xl0~ 2 


4.56xl0- 2 


1. 51xl0- 2 


9.27xl0" 3 


2.85xl0~ 3 


(f.ODXlU J 


(a \/ 1 n-2\ 
^4.04X1U J 


/ 1 c/i win - 2\ 
(1.04X1U j 


(y.oOXll) ) 


(z.ooxlU ) 


(2,0) - 


- (2,1) 


0.609 


0.522 


0.335 


0.253 


0.133 


(0.598) 


(0.511) 


(0.318) 


(0.250) 




(2,0) - 


- (3,0) 


4.45xl0~ 2 


2.86xl0~ 2 


9.69xlCT 3 


6.00xl0- 3 


1.90xl0" 3 


(4.50xl0- 2 ) 


(2.87xl0- 2 ) 


(9.74xl0- 3 ) 


(6.02xl0~ 3 ) 


(1.90xl0- 3 ) 


(3,0) - 


- (3,1) 


0.663 


0.582 


0.380 


0.308 


0.145 


(0.661) 


(0.577) 


(0.375) 


(0.299) 





Note. - Grid points are the same as in the tables of energy values. The values in 
parentheses are from 1DHF (Ruder et al. 1994). 



Table 9 

Oscillator strengths of bound-bound transition (m, 0) — > (m, 1) for carbon. 



TO 


P = 200 


P = 500 


p = 1000 





4.10xl0- 2 


1.30x10" 


-2 


5.81x10" 


-3 




(3.97x10-2) 


(1. 11x10" 


- 2 ) 


(5.71x10" 


- 3 ) 


1 


5.98x10-2 


2.25x10- 


-2 


1.28x10" 


-2 




(5.67x10-2) 


(2.31x10" 


") 


(1.35x10" 


- 2 ) 


2 


7.97x10-2 


3.92x10- 


-2 


2.47x10- 


-2 




(8.00x10-2) 


(3.80x10- 


- 2 ) 


(2.36x10- 


- 2 ) 


3 


0.114 


6.11x10" 


-2 


3.77x10- 


-2 




(0.109) 


(5.65x10- 


- 2 ) 


(3.62x10- 


- 2 ) 


4 


0.156 


8.63x10- 


-2 


5.72x10- 


-2 




(0.144) 


(7.94x10- 


2 ) 


(5.21x10" 


- 2 ) 


5 


0.197 


0.116 




7.83x10- 


-2 




(0.193) 


(0.113) 




(7.55x10- 


- 2 ) 



Note. - The values in parentheses are from 
1DHF (Miller 1990). The other electrons are in 
their ground state. 



